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Abstract 

The two-nucleon sector is near an infrared fixed point of QCD and as a result the S-wave scattering 
lengths are unnaturally large compared to the effective ranges and shape parameters. It is usually 
assumed that a lattice QCD simulation of the two-nucleon sector will require a lattice that is much 
larger than the scattering lengths in order to extract quantitative information. In this paper we 
point out that this does not have to be the case: lattice QCD simulations on much smaller lattices 
will produce rigorous results for nuclear physics. 
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I. INTRODUCTION 



One of the central goals of nuclear physics is to make rigorous predictions for both elastic and 
inelastic processes in multi-nucleon systems directly from QCD. The only presently-available 
technique to achieve this goal is lattice QCD, where space-time is discretized and QCD 
Green functions are evaluated in Euclidean space. Unfortunately, at present, the variety of 
processes that can be addressed with lattice QCD is quite limited. The currently-available 
computational power restricts not only the sizes of lattices that can be utilized, but also 
the lattice spacings and quark masses that can be simulated. Moreover, the Maiani- Testa 
theorem P| precludes determination of scattering amplitudes away from kinematic thresholds 
from Euclidean-space Green functions at infinite volume. However, by generalizing a result 
from non-relativistic quantum mechanics to quantum field theory, Liischer [3, Q realized 
that one can access 2 — 2 scattering amplitudes from lattice simulations performed at finite 
volume. Significant progress has been made using this finite-volume technique to determine 
the low-energy vrvr phase shifts directly from QCD, e.g. Ref. j^. However, only one lattice 
QCD calculation of the nucleon-nucleon (NN) scattering lengths has been attempted, 
and it was a quenched simulation with heavy pions ^. 

When contemplating computing nuclear observables with lattice QCD one naively as- 
sumes that the lattice must be much larger than the systems being simulated, so that the 
systems on the lattice resemble those at infinite-volume. This would mean, for instance, that 
when computing the rate for the simplest inelastic nuclear process, np —>■ dj, which near 
threshold involves radiative capture from the ^Sq channel, a lattice of size L ^ |a'''^^°^ | , |a*^^'^i^ | 
is required, where a^^^°^ and a'-^'^^-' are the ^5*0 and ^Si NN scattering lengths, respectively. 
Given that a^^°^ = —23.714 fm, such a calculation would have to await a future in which 
computational power is sufficient to handle volumes of this size. Fortunately, as we will see, 
this argument is not correct. 

There is a sizable separation of length scales in nuclear physics, due to the fact that 
nature has chosen to be very near an infrared fixed point of QCD p, |^ [l^. As a result, 
the scattering lengths in both S*— wave channels are unnaturally-large compared to all typ- 
ical strong-interaction length scales, including the range of the nuclear potential which is 
determined by the pion Compton wavelength. Perhaps counter-intuitively, in simulating 
two-nucleon processes, the relevant lengths scales are those of the nuclear potential and not 
the scattering lengths, and thus as long as the lattice is large compared to the inverse of the 
pion mass one can in principle "simply " determine matrix elements and scattering parame- 
ters. Furthermore, quantitative information about the two-nucleon sector can be extracted 
from simulations on even smaller lattices. However, the theoretical analysis that would be 
required for such an extraction is significantly more complex, and in this work we restrict 
ourselves to lattices that are much larger than the pion Compton wavelength ^. 



^ For a recent review of hadron-hadron interactions on the lattice, see Ref. [7|. 

^ There is a second technique that can be used to extract information about nuclear processes from small 
lattices. If simulations arc performed on lattices with quark masses that are somewhat different from their 
physical values, the scattering lengths will, most likely, no longer be unnaturally large [lli H^. By 
using the pionful effective field theory to calculate the quark mass dependence of the scattering parameters, 
the results of such simulations can be related to those at the physical values of the quark masses. The 
expansion parameters in the pionful theory are not exceptionally small, and so higher order calculations 
will need to be performed in order to have reliable extrapolations. We do not explore this option in this 
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Perhaps the motivation for a lattice calculation of the radiative-capture process np — >■ 
d'j is less than compelling as the cross-section and contributing multipoles at low-energies 
are well-measured. High-precision data recently collected 3 with HI7S agrees well with 
calculations jl^ [l^ in the pionless effective field theory 17, [l^, EFT(/), and also with 
the best modern potential models (see Ref. However, weak processes such as ud — >■ 

unp, which play a central role in the determination of solar-neutrino fluxes from the sun, 



depend upon two-body weak currents jl9 
uncertainties from reactor experiments 
tritium 



201 that presently are determined with significant 
_2ll | and are also determined from the /9-decay of 
2^ 1 with unknown systematic uncertainties. Recently they have been determined 
by the Sudbury Neutrino Observatory (SNO) from a fit to the neutrino fluxes ji^, but again 
with large uncertainty. As we move into an era of high-precision neutrino astronomy, and 
considering the potential impact this will have on our understanding of particle physics, it 
is imperative that we have a precision determination of these weak currents. Lattice QCD 
may be the only rigorous method with which to determine these capture rates with high 
precision unless a precise experimental determination is made 3 ■ 

In this work we explore the scattering states and bound states of the two-nucleon sector 
at finite lattice volumes. We first develop the finite-volume effective field theory relevant 
to the very-low energy interactions of two-nucleons, and we recover the exact eigenvalue 
equation as well as several approximate formulas due to Liischer jsi |^ . We find several new 
approximate formulas; we derive the leading finite-volume corrections to the bound-state 
energy and we find perturbative formulas for the lowest-lying energy levels for the case of a 
scattering length which is large compared to the lattice size. Armed with this technology, 
we consider simple unphysical limits of the scattering parameters to gain some intuition, 
and then explore the two-nucleon sector itself. 



II. THE PIONLESS THEORY OF NN INTERACTIONS IN A BOX 

An effective field theory (EFT) without pions, EFT(/), has been developed [1, H, [13, 113 to 
describe the very low-momentum interactions of two nucleons, with and without electroweak 
probes. EFT('^) exploits the sizable hierarchy between the 5— wave NN scattering lengths on 
the one hand, and the effective ranges, shape parameters and pion Compton wavelength on 
the other, and has been used successfully to perform relatively high-precision calculations, 
some at the ~ 1% level. Therefore, an EFT exists that can be used to rigorously determine 
the behavior of low-momentum nuclear observables in a finite volume, i.e. a lattice, and 
conversely can be used to extract infinite- volume limits of lattice simulations of two-nucleon 
observables. Furthermore, EFTj-j) has been used to successfully compute the properties 
of three-nucleon systems j25|, and therefore calculations of three-nucleon systems at finite 
volume should be possible as well ^. 

For NN scattering, the interaction between nucleons is described by a series of local 
operators with an increasing number of derivatives acting on the nucleon fields. The scat- 
tering amplitude can be computed in an elegant form using dimensional-regularization with 
power-divergence subtraction (PDS) [1,0]. In this scheme, the coefficients of the operators 



work. 

As a point of interest, it has been conjectured recently in Ref. |2g that QCD is very near the critical 
trajectory for a renormalization-group limit cycle in the three-nucleon sector. 
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in the Lagrange density have natural size, even for unnaturally-large scattering lengths. 
The scattering amplitude is identical to that of effective-r ang e theory, and that found by 
solving the Schrodinger equation with a pseudo-potential jl8|. An important feature that 
distinguishes EFT(^) from other constructions is that electroweak interactions can included 
systematically, in the same way they are included in chiral perturbation theory (^PT). 

In EFT(i^) (describing non-relativistic baryons ^ each of mass M) the exact two-body 
elastic scattering amplitude in the continuum arises from the diagrams shown in Fig. Q 
which can be resummed jsi |^ to give 




FIG. 1: The diagrams in EFT(^) which can he summed to give the scattering amplitude. The small 
solid circle denotes an insertion of the infinite tower of contact operators, X]C'2n(M) P^"" ■ 



where the C2n(/^) are the renormalization-scale dependent coefficients of operators with 2n 
derivatives acting on the nucleon fields (or equivalently with n time derivatives), fi is the 
dimensional-regularization scale, and D is the number of space-time dimensions. The loo 
integral Iq is linearly divergent, and when defined with the PDS subtraction scheme jT 
becomes 

ir'^ = {f^ + ^p) + OiD-4) , (2) 



where p = V ME is the momentum of each nucleon in the center-of-mass, and hence the 
scattering amphtude takes the usual form 

A=^ ^ . (3) 

M p cot 6 — ip 

This unitary expression describes NN scattering below the onset of the first inelastic thresh- 
old; that is, it is valid for |p| < a/ttv^? where m,r is the pion mass. The subtraction-scale 
dependence of the one-loop diagrams is exactly compensated by the corresponding depen- 
dence of the coefficients C2n(/^)- S is the energy-dependent 5— wave phase shift (we will only 
consider S"— wave scattering but this construction generalizes to all partial waves). It is clear 
from eq. (0) that p cot 6 is an analytic function of for momenta less than the cut-off of 
EFT(^), which is ~ m.,^/2. We may therefore adopt the effective-range expansion, 

11 °° 

PCoiS = + -ro / ' (4) 



Relativistic corrections can be included in perturbation theory 
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where a is the scattering length, tq is the effective range, and the other rj correspond to 
higher-order shape parameters ^. 

We are interested in the energy-eigenvalues of the NN system placed in a box with sides of 
length L with periodic boundary conditions. The scattering-state and bound-state energy- 
eigenvalues can be found by requiring the real part of the inverse scattering amplitude 
computed in the box to vanish. 



EC'2n(/i) P 



,2n 



Re{ir^>{L)) = , (5) 



where the infinite- volume integral in eq. (0) is replaced by a discrete sum over the momentum 
states allowed on the lattice, 

^oW = ^E^^ • (6) 

This discrete sum is also linearly divergent as the ultra-violet behavior of the theory is 
unchanged, and its value in the PDS scheme is found by adding and subtracting the corre- 
sponding infinite- volume integrals evaluated at E = 0. One of the infinite- volume integrals 
is evaluated with a momentum cut-off, |k| < A, that is equal to the mode cut-off introduced 
to regulate the discrete sum, while the other is evaluated with dimensional-regularization 
and PDS, to give 

^ ^ M V 7 I I 

and the limit A — > cxd is taken, assuming that the lattice spacing vanishes. For any realistic 
simulation there will be an upper bound on A, given by the edge of the first Brillouin 
zone i27il. 



A. The Eigenvalue Equation 

The energies of the low-lying energy levels of two-nucleons in a box with sides of length L 
with periodic boundary conditions 0, and with their center-of-mass at rest can 

now be determined in terms of pcot(5, from eq. © and eq. ((7j). Values of that solve ^ 

pcot.(rt^-is(gn , (8) 



^ We use the sign convention for the scattering length that is traditionally used in nuclear physics. This is 

opposite to the sign convention used by Liischer 
^ Liischer writes this expression as 



e 



2^So(k) ^ Zoo(l;g2)+Z^3/2^ 

Zoo(l; cp) - i7r3/2g 
where q — pi/(27r). The three-dimensional zeta- functions are 



where the formally divergent functions are defined via analytic continuation. 
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with 



1 



(9) 



IJI 



give the location of all of the energy-eigenstates in the box, including the bound states (with 
p2 < 0). The sum is over all three- vectors of integers j such that | j | < Aj and where the limit 
Aj — > oo is implicit (corresponding to the A — oo limit in eq. ((Tj)). In Fig. El we plot 8(77) 
vs. rj from eq. Q. While our derivation of eq. (jH)) is valid within the radius of convergence 
of EFT(^), that is for |p| < m.„/2, we expect that eq. (jSJ remains valid as long as the energy 
of the two-nucleon states are below the pion-production threshold, |p| < VrTvM 




FIG. 2: A plot o/S(r/) vs. rj from eq. The function has poles for r] > and does not have poles 
for rj < 0. 



B. Approximate Formulas 

There are two extreme limits that can be considered for the solution of eq. (jH)). First, there 
is the limit that Liischer considers in his work in which L ^ \a\. In this limit the solution 
of eq. (jHl) smoothly approaches the infinite- volume limit. The energy of the two lowest-lying 
continuum states in the Ai representation of the cubic group are 




(10) 



where the coefficients are Ci 



2.837297, Cs = +6.375183, and 



47r2 



12tan(5o r 



1 + c'itan(5o + C2tan^(5o + ... + 0{L~^) 



(11) 



ML2 



ML2 L 
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where c'^ = —0.061367, = —0.354156. In addition, in the hmit L ^ a, we have solved 
eq. (jHI) for the location of the bound state that exists for a > with an attractive interaction ^ 



7' 



1 + 11 I + 

jL l-27(pcot5)' 



(13) 



where (pcotS)' = pcot5 evaluated at = —7^. The quantity 7 is the solution of 



7 + pcot5\p2=_^2 = , (14) 

which yields the bound-state binding energy in the infinite- volume limit. 

In the limit where L <^ \{pcotS)~^\ (which is a useful limit to consider when systems 
have unnaturally- large scattering lengths), the solution of eq. (jH)) gives the energy of the 
lowest-lying state to be 

Eq = 'JJJji'^^ + C?2^PCOt(5o + ... ] , (15) 

where the coefficients are di = —0.095901, d2 = +0.0253716 and where p cot 60 is evaluated 
at an energy E = di. The energy of the next level is 

47r2 

El = -j^ [ d[ + 4 Lpcot^o + - ] , (16) 

where d[ = +0.472895, d'^ = +0.0790234 and where pcot^o is evaluated at an energy 

E = d[. The values of the df^ are determined by zeroes of the three-dimensional 

zeta-functions, and the expressions for Ei and Ei, excluding are valid for both-sign 
scattering lengths. 



C. A Toy Model : a = ±1 and = 

Let us consider an unphysical limit of NN scattering in which the NN potential has zero-range 
but a scattering length of \a\ = ±1 in the infinite-volume limit. Therefore, the scattering 
amplitude is pcot6 = — 1/a, as the effective range and all shape parameters vanish, = 0. 

The system with a = — 1 must result from an attractive interaction, but one not attractive 
enough to yield a bound state. One might imagine that the potential is extremely attractive 
and that the state with a = — 1 is the one near threshold with many other deep states 
present. However, the deep states will be at the cut-off of the theory, set by the range of 
the potential, and in the limit we are considering these are infinitely deep. The system with 
a = +1 could result from a repulsive interaction in which case there will be no bound state 
in this channel. However, a = +1 could also result from an attractive interaction that is 
attractive enough to give rise to a bound state near threshold. In the infinite-volume limit of 



^ The extension of the Chowla-Selberg formula to higher dimensions [s^l gives 

S{-x^) -2n^x + 67re-2'^=" + ... , (12) 

for large x, where the ellipses denote terms exponentially suppressed by factors of e^'^^^, or more. 
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FIG. 3: The two lowest-lying solutions to eq. ^ for a = ±1 fm and ri = 0. The vertical axis 
is , which is related to the energy by E = Q^-^jji, while the horizontal axis is logi^L. The left 
panel corresponds to a = —1 which can only arise from an attractive potential. The right panel 
corresponds to a = +1 which can arise from both an attractive and a repulsive potential. For a 
repulsive potential the lower solution is absent. The solid circles correspond to exact numerical 
solutions of eq. ^ . The curves that match the exact solution at small L result from eq. U5\) and 
eq. Ub}) . while the curves that match the exact solution at large L result from eq. U(J\) . eq. Ul\) . 
and eq. il.^) . 



this second scenario, one must recover the continuum of scattering states at positive energy 
and also the bound state. 

As L — »• the scattering lengths of opposite signs can be identified with each other (as this 
is equivalent to taking the |a| — >■ cxd limit of the model) and so the levels become degenerate. 
One is tempted to think that the levels that are degenerate in the L — >■ limit are the same 
as those that are degenerate in the L —>■ oo limit. However, this cannot be the case as the 
lowest state at L = in the system with a = +1 smoothly becomes the bound state at 
L = oo (as a function of L) and a bound state does not exist in the a = — 1 system. This 
can be seen clearly in Fig. |21 

It is apparent in both the L ^ |a| and L <^\a\ limits that the parameter arising in the 
asymptotic expansion of the energies of the levels is ~ 3|a|/L. This can be seen clearly in 
Fig. 121 where Liischer 's expressions m a/L break down at L ~ 3 fm while the expressions in 
eq. ()15p and eq. ()16|). which are an expansion in L/a, remain close to the exact solution out 
to L ~ 3 fm. 



D. Low-energy NN Scattering in the S-wave 

It is well known that the scattering lengths and effective ranges alone are sufficient to describe 
low-energy NN scattering data to quite high precision. This results in part from the fact 
that the shape-parameters are much smaller than one would naively guess. We therefore 
truncate the effective-range expansion in our numerical analysis, and the power-counting of 
EFT(^) dictates how the shape-parameters can be included in perturbation theory. At this 



8 



order, the energy-levels of two-nucleons in the ^5*0 channel whose center-of-mass is at rest in 
a periodic box of size L are found by solving 

' v^"^ + f r^n ^0 , (17) 



a(^5o) 2 nL \\2nJ 

for p^, which is related to the energy of the NN system via E = p^/M. In the ^5*0 channel 
the scattering length and effective range are 

^eso) ^ -23.714 fm , r^''^°'> = 2.734 fm . (18) 
Despite the fact that the NN interaction is attractive in this channel at long- and interme- 
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FIG. 4: The two lowest-lying energy- eigenstates of two nucleons in the Sq channel on a lattice of 
size L. The vertical axis is , where E = q'^jjj^, while the horizontal axis is the lattice size L. 
The solid circles correspond to the exact solution of eq. \n\ ) . The curves that are asymptotic to 
the exact solution at large L correspond to Liischer's relations in eq. U(J\) . eq. Ul\) . while the curves 
that are asymptotic to the exact solution at smaller values of L correspond to the expressions in 
eq. I115\) and eq. Illb]) . 



diate distances, there are no bound states of pp, np nor nn in the channel at the physical 
values of the quark masses or Aqcd ^■ 

One can see from Fig. 0] that Liischer's power-series expansion of the energy-levels 
converges slowly in the ^5*0 channel. For the ground state one needs a box of size 



One finds that in all likelihood there is no bound state in this channel for quark masses smaller than their 
physical values but it is possible that bound states exists for quark masses somewhat larger than their 
physical values This is an exciting possibility that can be explored with lattice QCD. 
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L> 80 fm ~ 3a^^^°^ while for the first excited state one needs a box of size L> 150 fm 
before these perturbative expansions converge to the exact result. Therefore, these asymp- 
totic expressions will not be of great utility to nuclear physicists in the near future. By 
contrast, the expressions we have derived in the L —>■ limit, eq. ()15|1 and eq. (fTT)|) . are 
applicable for boxes smaller than L< 50 fm. However, the crucial assumption, L ^ r,, 
begins to break down for volumes with L< 5 fm. In table H] we show the momenta of the 
first two states in the ^5*0 channel. The lowest-level can be described by EFT(^) on a lattice 
with L ~ 10 fm but a lattice with L> 15 fm is required in order for EFT(^) to describe 
the second state. This is not to say that we cannot use the location of all the states on 
a lattice with L ~ 10 fm; for momenta outside the range of validity of EFT(^) and below 
pion-production threshold, i.e. for ■m.„/2 < |p| < yJuiT^M, lattice results will have to be 
matched directly to p cot 5 in the ^5*0 channel as there is no effective-range expansion. 



TABLE I: Momenta of the lowest-lying levels of two-nucleons on the lattice. An asterisk denotes 
momenta outside the range of validity of the effective-range expansion, I p|™^^' = mT^jl ~ 70 MeV. 
The energy of the state is -E = |pp/Af, and the appearance of an "i" indicates a —ve energy. 





% |p 


(MeV) 


IpI (MeV) 


Lattice Size L (fm) 


1st 


2nd 


Deuteron 


1st 


1000 


0.1 i 


1.21 


45.5 i 


0.052 


100 


2.6 i 


10.52 


45.5 i 


1.76 


25 


13.8 i 


39.3 


45.8 i 


18.25 


15 


24.6 i 


67.0 


49.9 i 


44.61 


10 


39.0 i 


104.3 (*) 


61.3 i 


83.1 (*) 


5 


94.4 i (*) 


224.7 (*) 


116.5 i (*) 


206.5 (*) 



The '^S\ —^Di channel is somewhat complicated by the fact that the tensor interaction 
gives rise to mixing between the '^Si and the ^Di channels. On the lattice this means 
that different representations of the cubic group will mix due to the tensor component 
of the NN interaction. In particular, the Ai representation will mix with the E and T2 
representations [s^. However, the power-counting of EFT(^) dictates that we can ignore 
contributions from this mixing and the ^Di channel at the order to which we are working • 
Therefore, at this order, the energy levels of two nucleons in the ^Si channel can be found 
by solving 



1 




In the ^5*1 channel the effective-range parameters are 



+5.425 fm 







1.75 fm . 



(19) 



(20) 



One can see from Fig. that Liischer's power-series expansion of the energy-levels con- 
verges to the exact solution for the lowest-lying continuum state when L> 15 fm ~ Sa^^'^^^ 
However, our expression for the deuteron binding energy at finite-volume appears to work 
even at significantly smaller volumes. Unfortunately, the expressions that we have derived 
for small volumes do not converge well to the exact solution in this channel. This is be- 
cause the scattering length is only a factor of ~ 3 larger than the effective range and thus 
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FIG. 5: The two lowest-lying energy-eigenstates of two nucleons in the ^Si channel on a lattice of 
size L. The vertical axis is cp' , where E = q'^jjj^, while the horizontal axis is L. The solid circles 
correspond to exact solutions of eq. Uy\) . The curves that are asymptotic to the exact solution at 
large L correspond to Liischer's relation in eq. il(J\) and the relation for the bound-state energy, 
eq. m,'^). while the curves that are asymptotic to the exact solution at small L correspond to the 
expressions in eq. 1115]) and eq. ilh]) . 



the expansion parameter Lp cot S is not small enough over the entire range of L. Perhaps 
higher-order contributions will improve the agreement. The deuteron state can be described 
with EFT(^) for lattices with L> 10 fm, but the lowest-lying continuum state can be de- 
scribed by EFT(^) only for L> 15 fm. In table |l] we show the momenta of the deuteron 
bound state and the lowest-lying scattering state in the ^Si channel. 

Beyond the range of validity of EFT(^), the ^Si channel is significantly different than the 
^Sq channel. As mentioned above, in EFT(^) ^Si —^Di mixing is subleading in the expansion. 
However in the pionful theory, i.e. for < |p| < ^/m^rM, ^Si —^Di mixing appears at 

leading order in the EFT expansion js^l- Therefore, in this range of energies eq. (|19|) is not 
valid and one must solve a coupled system of integral equations. 



III. CONCLUSIONS 

Lattice QCD calculations of the scattering lengths and effective ranges in the two-nucleon 
sector would be a significant milestone toward rigorous calculations of nuclear properties 
and decays. Aside from providing essential information about the quark-mass dependence 
of nuclear physics, such calculations by themselves will not significantly improve our ability 
to compute nuclear properties, as the scattering amplitudes are already well-known exper- 
imentally. However, we would be in a position to compute electroweak matrix elements 
between two-nucleon states, and thereby provide information that will be vital to future 
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calculations of electroweak processes involving nuclei. Knowledge of these processes will 
directly affect analysis of data from present and future neutrino observatories. 

Exact solutions to Liischer's general formula for the energy-levels of the two-nucleon sys- 
tem on a lattice with periodic boundary conditions will allow for the extraction of scattering 
parameters from simulations with lattice volumes that are much smaller than naively esti- 
mated. It would appear that simulations on lattices with L> 15 fm will make it possible 
to extract both the scattering lengths and effective ranges in the two-nucleon sector in a 
straightforward way. This is contrary to the expectation that lattices with L ^ \a\ are 
required to determine the scattering parameters. The extraction of useful information from 
simulations with lattices L< 10 fm will require direct matching to pcot5 in the spin singlet 
channel, as the pionless theory and thus the effective-range expansion, will no longer be ap- 
propriate. In the spin-triplet channel there remains the additional challenge of formulating 
the finite-volume eigenvalue equations which account for mixing between S— and D— waves. 
These finite-volume calculations are a vital component of the technology required to make 
rigorous statements about nuclear processes directly from QCD. 
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